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' During the last decade, various modified theories of gravity have become very 

■ popular, f{R) gravity theory probably being the most studied one!^ Although most 

I models that are in agreement with observations are very close to general relativity, 

we have developed a much deeper understanding of the theory we wanted to mod- 
ify. During the year 2009 it was Hof ava-Lifshitz gravitjJ^ that excited the scientific 
I community with a new paper appearing on the subject area every other day. Hof ava 

proposed a power counting renormalizable theory for (3+l)-dimensional quantum 
gravity, which reduces to Einstein gravity with a non-vanishing cosmological con- 
stant in IR, but possesses improved UV behaviors. 

In this work, we explore the stability of the Einstein static universe in such 
I modified theories of gravity, for this study in general relativity see.^^ This can be 

motivated from various points of view. From a cosmological viewpoint it is the pos- 
sibility that the universe might have started out in an asymptotically Einstein static 
state, in the inflationary universe context.'^ On the other hand, the Einstein cosmos 
has always been of great interest in various gravitational theories. In general rela- 
tivity for instance, generalizations with non-constant pressure have been analyzed 
jjjjMl In tiie context of brane world models the Einstein static universe was investi- 
gated in,'2l while its generalization within Einstein-Cartan theory can be found in.^^^ 
Finally, in the context of loop quantum cosmology, we refer the reader to.'^^^ For 
the Einstein static universe in modified Gauss-Bonnet gravity, see .'^Finally, stabil- 
ity of the Einstein static universe in Hofava-Lifshitz gravity was analyzed in.'^HH 

By analyzing a simple background model and its perturbations one can study 
many properties of the modified theory in a rather explicit fashion. Let us give 
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one example to elucidate this point. When f{R) gravity became popular, it was 
believed that modifications of general relativity cannot stabilize solutions. However, 
as we showed explicitly in,^- this is not true and one can construct situations where 
the Einstein static universe, for instance, is stable with respect to a homogeneous 
perturbation. As we further showedj^^ this result in fact holds for all non-degenerate 
f{R) {f"{R) ^ etc) theories, see also. 1^ Generic stability results in f{R) gravity 
have been know since IDSS-^iSl 

The mathematics required to perform these investigations is straightforward and 
well understood, namely linear perturbation theory. The principal idea is to express 
all quantities in the form 

u\x) = u1g(x) +eu\{x) +e^u\{x) + . . . , (1) 

where the u\^J^x) describe a known exact solution of the field equations. It should 
be noted that there are situations where knowledge of u\^J^x) is not even necessary 
to solve to perturbed equations. In case of the Einstein static universe x would be 
the cosmological time t, u*(x) — {a{t) , p{t) , p{t)} , with a{t), p{t) and p{t) being 
the scale factor, the energy density and the pressure respectively. The background 
solution would be a static solution of the field equations which effectively reduce 
to algebraic equations that determine the background values. In general, we can 
assume the field equations to take the form 

F\u\ \7^u\ . . . , Va^VpV^Vsu') = 0. (2) 

Here we restrict ourselves to theories which contain up to four derivatives. We 
also assume that the number of field equations matches the numbers of unknown 
functions, we exclude over-determined and under-determined systems. Notice that 
the field equations sourced by a perfect fiuid are under-determined as long as no 
equation of state is assumed. Therefore, we consider a perfect fluid (no anisotropic 
stresses) with a linear barotropic equation of state, Pbg(pbg) = wpbg- Now, ii p{p) = 
wp is also assumed to hold at perturbative level then the sound speed of adiabatic 
pressure perturbations is given as — w. Note that it would be interesting to 
investigate the situation when this assumption is dropped, ^ w. 

The approach described covers usual General Relativity and thus the entire field 
of cosmological perturbation theory, f{R) gravity, modified Gauss-Bonnet gravity, 
or f{G) gi'avity, and also Hofava-Lifshitz gravity. Note that the above treatment 
is sufficiently general to also include brane world models which pose additional 
technical challenges .'^Slsl] 

In order to investigate the behavior of the perturbation, one now inserts Eq. ([IJ 
into the field equations ([2]) and linearizes the equations with respect to e (Taylor 
expansion about e = 0). In full generality, the field equations then become 

F'«g, V„<g, . . . , V„V^V^V5<g) 

+ < a^j bg + • • ■ + ^ ( av.v,v,v..- ) ,g V" + - «• 
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Those equations are now linear in the perturbed variables. In Gauss-Bonnet gravity, 
for example, this equation takes the simple form 

24^3^2^(1 + wff"{0) a'l"{t) + 2a'i{t) ~ Kpbg(l + w){l + 3u)) ai{t) = 0. (4) 

In this equation ai is the perturbed scale factor and / is the function which deter- 
mines the modifications of the Gauss-Bonnet term. 

In general one finds a set of linear, coupled differential equations. Equations 
of this type can always be solved analytically. We can therefore conclude that the 
additional degrees of freedom in any modified gravity model lead to enhanced regions 
of stability in the parameter space. 
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